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1. Lotka-Volterra competition and species distribution 

Competitive interactions occur when entities in a system grow by consuming com- 
mon finite resources. They are ubiquitous in many fields of science: examples in- 
clude biological species competing for food (MacArthur & Levins 1967; Roughgar- 
den 1979; Case 1981), mode competition in nonlinear optical systems (Benkert & 
Anderson 1991), or alternative technologies competing for a market (Pistorius & 
Utterback 1997). An early, simple, but powerful model for competitive interactions 
is the Lotka-Volterra (LV) set of competition equations (Volterra 1926; Lotka 1932): 




N i = r i N i \l- — yG ij N j \, i = l,...,m. (1.1) 



where m is the number of species, Ni the population of species i, its maximum 
growth rate, Ki its carrying capacity, and Gy is the matrix characterizing the 
interaction among species i and j, more specifically the decreasing on the growth 
rate of species i by the presence of j. Competitive interactions are characterized 
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by Gij > 0, the situation to be considered here, whereas negative interactions may 
model situations of mutualism, predation or symbiosis. 

In classical ecological niche theory, species are associated to points in an abstract 
niche space. Coordinates in this space represent relevant phenotypic characteristics, 
for example size of individuals in a species, or the size of preferred prey, such that 
intensity of competition is larger if species are closer in this space. We assume for 
simplicity this space to be one-dimensional (multi-dimensional generalizations are 
straightforward, as briefly mentioned later). If niche locations can be considered to 
be a continuum, we can write Eq. (1.1) as: 



where now ip(u, t) is the population density at niche location u. The integral extends 
over the full niche space, which could be finite or infinity. For most purposes, Eqs. 
(1.1) and (1.2) can be considered as equivalent, since the second is obtained from 
the first in the limit of many close interacting species, and (1.1) can be recovered 
from (1.2) for a discrete distribution of species: 



with = G(ui,Uj), Ti = r{ui) and = K{ui). 

It is widely believed that (1.1) or (1.2) predict a competitive exclusion leading 
to a limiting similarity situation (Abrams 1983), in which a pair of species too close 
in niche space can not coexist, and one of them would become extinct. However it is 
known that the model allows for continuous coexistence of species in some situations 
(Roughgardcn 1979), and refinements on the conditions for this coexistence have 
been developed, with emphasis on the effect of the shape of the carrying capacity 
function K(u) (Meszena et al. 2006; Szabo & Meszena 2006). In this context, a 
particulary surprising result was the finding by Scheffer & van Nes (2006) of a 
situation -for uniform carrying capacity- which was neither of full coexistence nor 
of full exclusion, but of clusters or lumps of tightly packed species which did not 
exclude each other, but were well separated from other clusters so that there was 
a type of limiting similarity leading to a minimum intercluster distance. Clustering 
of individuals or entities under competitive interactions of the LV type had been 
already observed in other contexts (Fuentcs et al. 2003; Hernandez-Garcfa & Lopez 
2004, 2005; Ramos et al. 2008), where the mechanism was the diffusive broadening 
of an otherwise zero-width species or entity. In contrast, the lumps in Scheffer & 
van Nes (2006) appeared even in the absence of diffusion in niche space, which is 
the situation also considered here. 

The importance of the functional form of the interaction kernel Gij in (1.1) or 
G(u, v) in (1.2) was stressed by Pigolotti et al. (2007) for the case of uniform carry- 
ing capacity and interactions depending only on differences of niche positions, and 
found to be relevant in an evolutionary context by Leimar et al. (2008). For that 
case the positive-definite character of the Fourier transform of G(u, v) = G(u — v) 
is a condition implying the absence of limiting similarity. Species clustering was 
reported, but for interaction functions rather different from the Gaussian used in 
Scheffer & van Nes (2006). In fact, for the Gaussian interaction case most results are 
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extremely sensitive to details such as the implementation of the boundary conditions 
or weak ecological second order effects (Pigolotti et al. 2008). Thus, a clarification 
of the mechanisms leading to species clustering in LV models would be desirable. 
In addition, the results in Pigolotti et al. (2007, 2008) were obtained under the 
unrealistic assumption of homogeneity in niche space whereas the inhomogeneities 
in the carrying capacity are known to play relevant roles (Szabo & Meszena 2006). 
For simplicity we restrict our description to the standard situation in which com- 
petition is stronger among species closer in niche space. It is worth mentioning 
the existence of studies of LV systems where non-local interactions are considered 
(Doebeli & Dieckmann 2000). That situation can also be described by the general 
formalism used here of an integral kernel function, and our general results therefore 
also apply to the situation with non-local interactions. 

In this Paper we analyse some mathematical properties of the LV model (1.1) 
or (1.2). In Sect. 2 we show that the positive-definiteness of the kernel G remains a 
determining condition for stable coexistence even for non constant K(u). In Sect. 
3, we discuss the mechanism producing lumped species distributions and explicitly 
give an analytic expression for a particular interaction kernel. In the Appendix we 
show that, in contrast with the earliest characterizations of the interaction kernel G 
(MacArthur & Levins 1967; Roughgarden 1979), both positive- and non-positive- 
definite kernels can arise from more detailed ecological models which consider the 
dynamics of the consumed resource. We use periodic boundary conditions in our 
numerical simulations. We expect the effects of this simplifying but unrealistic as- 
sumption to be unimportant at least when a non-constant carrying capacity limits 
the presence of species to a limited region of niche space. 

2. The stability of close coexistence 

A simplifying assumption for the study of the LV model is that of homogeneity in 
niche space. In this case, the carrying capacity and growth rate are constants, K 
and ro, and the interaction kernel depends only on differences of niche positions 
G(u, v) — G(\u—v\). Niche space could be infinite, but in the case in which it is finite, 
homogeneity can only be achieved under periodic boundary conditions. Under these 
restrictions it is easy to see that a steady solution to (1.2) which is homogeneous 
and everywhere non-vanishing always exists: V'o = Kq/Gq, where Go = / duG(u). 
This solution represents coexistence of all possible species without a limit to their 
similarity. Its stability against small perturbations can be analysed by linearization 
of the equation resulting from substitution of ip(u, t) = ipo + 6ip(u, t) into (1.2). The 
solution for the Fourier transform of the deviation from the homogeneous state, 
6tp q (t), is 

= S^ q (0)e x ^ .with \ q = -r ^ . (2.1) 

Go 

where G q is the Fourier transform of G{u). Thus, the homogeneous solution tpo 
is stable if G q is positive Vg, while a instability leading to pattern formation oc- 
curs when G q may take negative values (Pigolotti et al. 2007; Fuentes et al. 2004: 
Hernandez-Garcfa & Lopez 2004; Lopez & Hernandez-Garci'a 2004). We note that 
many steady solutions to Eq. (1.2) exist besides tpo (in particular, solutions of the 
form (1.3)). This is so because dynamics preserves ip(u) = at all places where 
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there is no initial population. Notice also that ip is the only strictly positive solu- 
tion. Among this multiplicity of solutions the ones that will be more relevant are 
the ones which are stable under perturbations or small immigration (Pigolotti et al. 
2007). 

An interesting class of functions to be used as kernels and carrying capacities is 
the family {g?} given by 

«£(«) = expHuH"), (2.2) 

which is parameterized by the value of p. The widely used Gaussian kernel is ob- 
tained for p = 2. When p < 2 the functions are more peaked around u = (the 
case p — 1 is an exponential) and for p > 2 they become more box- like (g^(u) is 
the flat box with value 1 in the interval [— a, a] and zero outside). The width of 
the kernel a gives the competition range in niche space. We have positivity of the 
Fourier transform if p < 2. This implies that the homogeneous solution is stable 
under evolution with uniform K and kernel G of the form (2.2) if p < 2. When 
p > 2, the homogeneous solution is unstable and the system approaches delta comb 
solutions of the type (1.3), with a spacing approximately 1.4<r (Pigolotti et al. 2007) 
which represent limiting similarity situations. 

We now generalize the above stability analysis to the more realistic case in 
which there is no homogeneity in niche space. First we consider the simpler case 
of a symmetric kernel G(u,v) = G(v,u), which in particular includes the previous 
case of kernels depending only on species distance: G(u, v) = G(\u — v\). Note that 
in this symmetric case one can write Eq. (1.2) in potential form: 

flrtKM) = -r(«) K{u) mu) t (2.3) 

with the functional potential given by: 

V[ip] = - J K(u)ip(u, t)du +\J J G K v)ip(u, t)tp(v, t) du dv. (2.4) 

Stationary solutions of Eq. (1.2) are those for which the r.h.s of Eq. (2.3) equals 
0. This has many possible solutions. We define the natural stationary solution, 
ip N (u), as the one which is positive and non vanishing for all u, so that 



SV\ 



= 0, (2.5) 



that is, the one satisfying: 

J G(u, v)ip N (v)dv = K(u) , V u (2.6) 

The solution tp N (u) can be considered the non-homogeneous generalization of ipo 
introduced in the homogeneous case. In the particular case in which G(u, v) — 
G(u—v) the natural solution can be explicitly written in terms of Fourier transforms 
of the competition kernel and the carrying capacity, either in an infinite system or 
in a finite one with periodic boundary conditions: 

V? = ^- (2.7) 

G q 
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Figure 1. Long-time solutions of (1.2) for different kernels and carrying capacities. Left: 
G = g\, K — sech(u/a), with a = 0.1. The natural steady solution (ip N = a~ 1 secb 3 (u/a)), 
which is positive and non-vanishing everywhere, is reached at long times. Center: G = go.ii 
K — flo.i- Under this non-positive-definite competition kernel, the solution shown is still 
evolving and approaches a singular delta comb of the type (1.3) at long times. Right: 
G — go.'ii K — go, i. A positive natural solution does not exist and the system approaches 
a single hump solution which vanishes in part of niche space. 



This requires that these Fourier transforms and their inverses exist and lead to 
positive populations densities. When this happens, a continuum species coexistence 
is obtained, and its existence is generally robust against small changes in G or K. 
We show later that it is also an attractor of the dynamics when G q satisfy positivity 
requirements (p < 2, for the family in (2.2), being p = 2 the marginal case). For a 
uniform carrying capacity, the natural solution (9) always exists and is uniform in 
phenotype space tp N (u) = tpo. But the natural solution may lose positivity or even 
cease to exist depending on the properties of G and K. For example, when both 
G(u) and K(u) are of the form (2.2) with p = 2, the inverse Fourier transform of 
(2.7) exists when the carrying capacity has a value of a larger than the kernel G, 
but not in the opposite case. 

Figure 1 shows stationary solutions attained at long times by the dynamics 
in (1.2) illustrating the situations described above, starting from a smooth non- 
vanishing initial condition. In the first case we choose a kernel and carrying capacity 
functions (G(u) = g],{u), K{u) = scch(u/a)) such that the natural solution exists 
and is positive everywhere. Thus it is stable, and it is the steady state attained at 
long times. In fact it can be analytically calculated: 

4> N {u) = a-^ech^u/a) . (2.8) 

In the second case the non-positiveness of the kernel used (with a carrying 
capacity of the type Eq.(2.2)) breaks down the initial configuration into lumps, 
which at long times approach zero-width delta functions with forbidden zones in 
between. In the third case, despite G q being positive, a positive natural solution 
does not exist. Several outcomes are possible but for the kernel and capacity used, 
the system approaches a single hump solution which vanishes in part of the niche 
space. 

More in general, but still in the symmetric G case G(u,v) = G(v,u), writing 
the LV model in potential form (Eq. (2.3)) is of great use since one can show 
that, provided r(u) and K{u) are positive, dV/dt < 0. This implies that V is 
a Lyapunov potential and dynamics proceeds towards its absolute minimum, or if 
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ip{u, t = 0) = for some u, towards the minimum of V under such constrain. Notice 
that, since the potential is a quadratic form, %jj N is a global attractor (starting 
from non-vanishing initial conditions) when the competition kernel is a positive 
definite quadratic form, which means that J J f(u)G(u,v)f(v)dudv > 0, V/ (or 
J2ij x-iGijXj > 0, y{xi\ in the discrete case). This generalizes the previous stability 

condition on the Fourier transform G q > to niche inhomogeneous cases, and shows 
that the stability result was global indeed. In a multi-dimensional niche space the 
same analysis shows that the positive-definiteness of the quadratic form remains the 
condition for the global stability of the natural solution. In any case, the important 
consequence is that the stability of the natural solution depends uniquely on the 
competition kernel and not on the carrying capacity (provided the relation kernel- 
capacity is such that the natural solution exists and is positive). In particular, for 
competition kernels of the form (2.2), ip N is always (if existing and positive) a 
globally stable solution of the dynamics for p < 2, and unstable otherwise. 

The crucial difference in the case of a non-symmetric competition kernel is that 
there is no obvious Lyapunov potential for the system. This implies that there are 
no available global stability results. However, local stability can be investigated. Let 
us consider a small perturbation of the positive natural solution ip N (u) + 5ip(u,t). 
To linear order, the perturbation evolves as: 

dSip(u,t) ,ip N (u) 



dt 



= -r(«)^ J G(u, v) 6^(v, t) dv. (2.9) 



We now consider the functional H(5ip) = J du (A(u)i ; T(u)/r(u)'i/' JV (w)) (Sip) 2 , 
where A(u) is a positive function so that H > and H (0) = 0. Let us compute its 
time derivative: 

^ = -2 J 6tp(u) A(u) G{u,v) 5tp(v)dudv. (2.10) 

If for some choice of A(u) one has that A(u)G(u, v) is positive definite, then 
dHj dt < and Sip = will be approached. This shows that ij) N is linearly stable in 
such case. We note that the case in which G(u, v) itself is positive-definite trivially 
guaranties the positivity of A(u)G(u, v), with a constant A. Thus, even in this more 
general nonsymmetric case, it is the character of the interaction kernel G, and not 
of the carrying capacity (provided it is such that the natural solution exists and is 
positive), which determines the stability of the natural solution. 



3. Lumped species distributions 

Schcffer & van Nes (2006) found transient but long-lived solutions of Eq. (1.1) con- 
sisting of periodically spaced lumps containing many close species. They used a 
Gaussian interaction kernel which turned out to introduce an excessive sensitivity 
of the results to the numerical implementation of the model and boundary condi- 
tions (Pigolotti et al. 2008). They found however similar solutions as steady con- 
figurations when adding an extra predation term acting effectively only on species 
with high population. This can be thought as an extra intraspeciGc competition 
since it decreases the growth of species with many individuals. Exploiting this idea, 
Pigolotti et al. (2007) checked the effect of using in (1.2) a kernel of the type (2.2) 
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but with an enhanced interaction at u = 0, i.e. enhanced intraspecific competition. 
In particular, they used a constant carrying capacity K(u) = Kq and a flat box 
kernel with an added delta function at the origin (see Fig. 2), 



G(u)=g™(u)+a5(u) 



(3.1) 



Lumped patterns were obtained numerically for a = 1. 

Because the dynamics of (1.2) usually involves very long transients, it is inter- 
esting to calculate analytically the steady lumped solution in the simple case of a 
kernel (3.1) and uniform carrying capacity Kq (in the infinite line). 
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Figure 2. The kernel in Eq. (3.1) (left), and the analytic steady solution given by (3.5) 
and (3.8-3.9) for a = K = 1, a = 0.8, L = 0.3 and d = 1 (right). 



We begin with the steady state condition 

J G(u,v)ijj(v)dv = K{u) , (3.2) 
valid at u such that ip(u) ^ 0, that particularized to (3.1) and constant K reads: 

aip(u) + / dvip(v) = K . (3.3) 

J u — a 

This is transformed into a differential-difference equation after differentiation with 
respect to u: 

aip'{u) + ip(u + a) - ip{u - a) = , where tp(u) ^ . (3.4) 

This steady equation has many solutions, including the natural one ipo = Ko/(a + 
2a) which is non- vanishing everywhere, or delta combs such as (1.3). We search for 
solutions of the type in Fig. 2, i.e. periodic arrays of lumps, of period d, each one 
having a symmetric hump shape f(u) of width 2L (i.e. f(u) — if u ^ [—L,L]): 

oo 

lK«)= E /(z-™0 ( 3 - 5 ) 

n— — oo 

We are assuming that the lumps do not overlap, so that d > 2L. We also note that 
if <T + 2L < d there is no interaction between different lumps, so that for u £ [-L, L] 
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Figure 3. The solution A (positive branch) of Eq. (3.7), giving the inverse width of species 
lumps. The width is finite for d — a < a, which is favored by larger enhanced intraspecific 
competition a. 

Eq. (3.4) reduces to f'(u) = and there is no lump solution. Moreover, analysis is 
much simplified if each of the lumps interacts only with its neighbors (a + 2L < 2d) . 
Thus we restrict to d < a + 2L < 2d, for which (3.4) with (3.5) and u € {—L, L) 
becomes: 

af'{u) + f(u + a - d) - /(it - a + d) = (3.6) 

The general solution of this linear equation is obtained as a sum of exponentials 
cxp(Au), with 

aX = sinh(A(d- a)) . (3.7) 

A = is always a solution, and if d — a < a there are two additional solutions 
±A, plotted in Fig. 3. For d — a > a the only solution is the constant one, but in 
the opposite case (the situation favored by enhanced intraspecific competition a) 
the solution is a linear combination of three exponentials. Two of the constants of 
the combination are determined from f(L) = f(—L) = 0. The third one, which 
gives the overall normalization, can be obtained by returning back to the original 
equation (3.3). Finally we get 

><«> - 'O--^)""'-^ 

= elsewhere (3.8) 

with 

A = ^ (3 Q) 

a(l-sech(AL)) + f (AL - tanh(AL)) ' y ' ' 

and the value of A which is plotted in Fig. 3. Figure 2 shows the analytic solution 
(3.5) with (3.8)-(3.9). We have not studied the stability of this configuration. But 
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the numerical results in Pigolotti et al. (2007) indicate that it is stable for some 
values of L and d. 



2.0^ 

1.5: 









-0.4 -0.2 



Figure 4. The kernel G = go. 2 + 0.8g 0.02 (left)) and the steady solution obtained 
numerically from it at long times with constant Ko = 1 (right). 

We finally stress that the appearance of the lumped solution is not a consequence 
of the singularity of the delta function in the kernel. In fact, any kernel sufficiently 
peaked at the origin will favor the coexistence of close species. If the behavior at 
larger distances of the kernel makes it not positive-definite, then full coexistence will 
be unstable and the natural solution will split into disjoint lumps. An example of the 
final steady state in this situation is shown in Fig. 4, with a kernel G = g$ 2 +0.8(?o 02 
which has the properties just described and contains no delta singularity. 

C.L. and E.H-G. acknowledge support from project FISICOS (FIS2007-60327) of MEC 
and FEDER and NEST-Complexity project PATRES (043268). K.H.A was supported by 
the Danish Research Council, grant no. 272-07-0485 



Appendix A. Models leading to LV competitive interactions 

We have seen that the character of the interaction kernel G is of major importance 
to determine the qualitative outcome of LV competition. In the original formula- 
tion of the niche model, however, only positive definite kernels were allowed. The 
reason is that competition kernels were derived in terms of utilization functions 
Ui(x), describing how consumer i uses resource at niche location x (assumed to be 
continuous) (MacArthur & Levins 1967; Roughgarden 1979): 



J Ui{x)uj(x) dx 
J uf{x) dx 



(Al) 



When the resource is directly related to space, (A 1) can be justified by considering 
the probability that consumer i meets consumer j (Roughgarden 1979). It is easy 
to see that Gij obtained from (A 1) is positive definite. We show in the following, 
however, that relation (A 1) is by no means general, and that a greater variety of 
kernels -positive or non-positive definite, so that the natural solution representing 
coexistence can be either stable or unstable- could be obtained from equations in 
which resources are explicitly modelled. Related calculations could be found, for 
example, in Schoener (1974). 

We consider a set of predators (or consumers), with populations Ni, i = 1,2, ...m, 
competing for different types of prey populations or resources, R a , a = 1,2, ...n, 
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the later growing in a logistic way with growth rate (3 a and carrying capacity Q a 
in the absence of predators. Particular equations modelling this are 

R a = -RaY^^i + ^Rc^-^j (A 2) 

Ni = N l J2S ia R a -d 2 N l (A3) 

a 

di is the death rate of species i. The interaction coefficients are a a i, the depiction 
rate of resource a produced by species i, and the sensitivity Si a , giving the growth 
rate of i thanks to resource a. Lotka-Volterra type dynamics arises when the time 
scale for resource evolution is much faster than that of the consumers (i.e. Si a 
and di — > oo, but with their ratio finite). In this case, adiabatic elimination of the 
resource can be done (R a « 0, so that each prey is at each instant at the equilibrium 
determined by their consumers), giving 

Ra^Qa(l-^-^2a ai N^j . (A 4) 

for the non- vanishing resources. The impact matrix, D ai , describing the depletion of 
resource a by species i (Meszena et al. 2006), is D a i = Q a a a i/ (3 a , which substituted 
into the consumers equation leads to: 

iVi=JVi ^ , X r " V j ' ( A5 ) 

where r; = J2 a Si a Q a is the maximum growth rate and Cij — ^2a SiaD a j. Thus, 
the result is an effective interaction among the predators which is of Lotka-Volterra 
type. It is customary to write (A 5) in terms of the carrying capacity Ki, defined 
as the equilibrium population Ni attained in the absence of the other competitors, 
i.e. Ki — ri/Cu. In terms of it, Eq. (A 5) becomes identical to (1.1), with 

_ Cij _ J2a SjgDaj . . . 

°« ~ ~n7. ~ v <? D ' [ ' 

Having a continuum R(x) of resources instead of a discrete set R a does not intro- 
duce important difficulties. Simply one should replace sums by integrals, replacing 
the coefficients of Eq. (A 5) by: 

n = J ' S l {x)Q{x)dx , (A 7) 

Cij = j Si(x)Dj(x)dx , (A 8) 

J S l (x)D :j (x)dx 
ij ~ J Si(x)Di(x)dx ' 1 j 

One can also consider a continuum of species, labelled by their phenotypes u, 
so that Eq. (1.1) is replaced by Eq. (1.2) with K(u) = r(u)/C(u,u), G(u,v) — 
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C(u, v)/C(u, u), and r(u) and C(u, v) given by obvious generalizations of (A 7) and 
(A8)._ 

It is clear that the presence in the kernel of two different functions (compare 
with the most restrictive expression (A 1)) gives enough freedom to obtain a variety 
of kernel behaviors under different circumstances. A particularly popular choice is 
to assume that impact and sensitivity are proportional: Si a = eD a i, with a constant 
efficiency e. In the continuum resource case the functions can be written in terms of 
a single utilization function Ui(x) as Di{x) = Ui(x) and 5, (x) = eui(x), leading to 
the classical expression (A 1). Slightly more general cases arise when the efficiency 
depends only on the resource, e = e a , or on the consumer e = ej, or when dependence 
on the two types of species factorizes, e = ViW a . In all these cases (if the efficiency 
is positive) one is lead to a kernel Gij which is positive definite. In more general 
cases, one can have a kernel leading to instability of the coexistence state. 

We conclude with two instances of ecological interactions in (A3) which allow 
to tune the stability. First, a homogeneous discrete and infinite niche space in which 
all resources have the same internal dynamics Q a = Q, (3 a = (3, Ma, as well as the 
consumers: rf, = d, Mi. The interactions are taken to be 

Sia = gdi, a (A 10) 

0-a,j = ^® a i = a $a,j + K^oj-l + 5 a ,j + l)- 

This models a situation in which the consumer k grows only by consuming its 
optimal resource Rk, whereas it depletes also the neighboring resources, Rk+i and 
Rk-i- We have n = Qg, K { = (3/a, Qj = (Qg/(3) (aS itj + b(S itj - 1 + S itj+1 )), and 
Gij = Sij + (b/a) (Sij-i + Sij + i) so that equation (1.1) is now 



Ni - QgNi 



l-^iaNi + biNi+i+Ni-i)) 



(All) 



The natural solution, i.e. the one in which all species have positive non zero pop- 
ulation, is Ni — (3/(a + 26), Mi. Its linear stability can be studied by linearization, 
Ni{t) = N t + 8Ni(t) and substitution of the ansatz SNi w e x ^e iql (here i = \/ z T). 
We find A g = —{Qg/P) (a + 26cosg), q G [— tt,it]. X q are the eigenvalues of — Cij, 
and stability of Ni requires all these eigenvalues to be negative, i.e., Cij to be pos- 
itive definite. When a > 26, then A g < Vq, and the natural coexistence solution 
is globally stable (see results in Sect. 2). It is unstable otherwise. In this example, 
there is no well-defined single utilization function and the positivity properties of 
the interaction kernel and thus the stability of the natural solution can be changed 
by varying the parameters. 

As a second example, with nonconstant carrying capacity, we consider a contin- 
uous distribution of resources and species on the line, and we take 

Q(x) = Qg(x), (3{x) = (3g{x), 
S u (x) = s5(u-x), a u (x) = f(u-x). 

which implies that the consumers of phenotype u grow only from the resource at 
location x — u, but they deplete a wider range characterized by /. This leads (in 
the continuous formalism) to r{u) = sQg(u) , C(u,v) = sQf(u — v)/(3, K(u) = 
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(3g(u)/f(0) and G(u,v) = f(u — v)/f(0). In this example, by choosing the func- 
tions g and /, we can impose any desired combinations of carrying capacity and 
interaction kernel. Gaussianity or positive definiteness are particular cases, no more 
natural in this generalization than alternative choices leading to non-positivcness, 
instability, and thus exclusion zones between clumps of species. 
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